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Abstract 



We investigate the nucleation of superconductivity in a uniform perpendicular 
magnetic field H in aluminum microsquares containing a few (2 and 4) submi- 
cron holes (antidots). The normal/superconducting phase boundary Tc{H) 
of these structures shows a quite different behavior in low and high fields. 
In the low magnetic field regime fluxoid quantization around each antidot 
leads to oscillations in Tc{H), expected from the specific sample geometry, 
and reminiscent of the network behavior. In high magnetic fields, the Tc{H) 
boundaries of the perforated and a reference non-perforated microsquare re- 
veal cusps at the same values of $/^'o (where ^ is the applied flux threading 
the total square area and $o is the superconducting flux quantum), while 
the background on Tc{H) becomes quasi-linear, indicating that a giant vor- 
tex state is established. The influence of the actual geometries on Tc{H) is 
analyzed in the framework of the linearized Ginzburg-Landau theory. 
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I. INTRODUCTION 



Recent experiments on mesoscopic superconducting aluminum structures with sizes 
smaller than the temperature dependent coherence length ^{T) and penetration depth A(T) 
have shown the influence of the sample topology on the superconducting critical parameters, 
like the normal/superconducting phase boundary Tc{H) (i.e. the critical temperature Tc in 
the presence of an applied mag netic field H)^. Many different topologies have been studied 
experimentally and theoretically. 

First of all, there are structures made of quasi-one-dimensional strips, which can be 
further classified into sing le loopM (where Tc{H) shows the well-known periodic Little- 
Parks oscillations!), multiloop structureJl^ (double loop, yin-yang, lasso, 2x2-cell, etc.) 
and large infinite network The theories used to calculate the Tc{H) for these structures 
are based on the linearized Ginzburg-Landau theory, using either the de Gennes- Alexander 
formalism§ or the London limilJi§. The fiuxoid quantization constraint (i.e. the require- 
ment that the order parameter is uniquely defined after integration of the phase gradient 
along a closed superconducting contour) in these 'multiply connected' structures gives rise 
to the oscillatory shape of the phase boundaries Tc{H), superimposed usually on a parabolic 
background. 

Secondly, surface superconductivity effects in (circular or square shaped) single dot JiUiHii 
and antidotdll^ (i.e. one antidot in a plain film) have been studied intensively. In 
these structures, Tc{H) consists of oscillations, which are pseudoperiodic. The appear- 
ance of thegiant vortex state, where superconductivity is nucleated only near the sample 
boundarytJii^MEHnHHl, is due to the quantization of the phase winding number L of the super- 
conducting order parameter \1/ = e"*'^'^ (this is equivalent to fiuxoid quantization). For 
cylindrically symmetric structures, one refers to L as the angular (or orbital) momentum 
quantum numberi. For the states L > 1, the dot (or antidot) area is threaded by multiples 
L of the superconducting flux quantum $o = h/2e. This "surface superconductivity" gives 
rise to a quasi-linear critical field Hcs versus temperature T, which we will further compare 
with the experimental Tc{H). 

The experimental studies of the antidot structure are usually carried out on samples 
with regular lattices of antidots. It was shown by Bezryadin et al. that, at sufficiently high 
magnetic fields, the antidots behave independent lylil. In these systems, the antidots create 
efficient pinning centers for the flux line latticeii. 

We will present the measured phase boundaries Tc{H) of three different topologies, which 
are shown in Fig. |I[ The three structures studied are a filled microsquare, and two squares 
with 2 and 4 square antidots respectively. Similar structures were studied in Refs.i'i, where 
the 4-antidot structure was proposed as a basic cell for a memory based on flux logic. In 
those papers, different stable vortex configurations were detected at low magnetic fields. 

The goal of the present report is to study the influence of the antidots inside a microsquare 
on the crossover from the "network" behavior at low fields to a giant vortex stateillEi at 
higher fields, and whether eventually the two configurations (vortices pinned by the antidots 
and the giant vortex state) can coexist. We will mainly focus on the high magnetic field 
regime. In a structure with one single antidot (i.e. a loop, as in Ref.0) as well, the presence 
of a giant vortex state can be anticipated at sufficiently high magnetic fields. In such a 
system, however, the phase winding number L will be identical for every contour encircling 



2 



the antidot. The development of the giant vortex state is accompanied by a transition from 
a parabohc background in Tc{H) to a quasi-hnear Tc{H) behavior. The crossover field is 
strongly dependent on the size and the aspect ratio of the loop and will be the subject of a 
future paper. For the loop studied in Ref.0, the magnetic field was clearly not sufficiently 
high to reach this transition regime. 

The advantage of a structure with more than one antidot is the property that each 
antidot can in principle contain a different number L of flux quanta $o- Simultaneously, a 
quantum number L is attributed to the outside square. The observed cusps in Tc{H) can 
then be related to the switching of either the quantum state of an antidot, or of the whole 
square. 

For the 4-antidot structure, a 'collective' or network behavior can be expected at low 
magnetic fieldsl, while a 'single object' regimeEl can be reached at higher fields, where at 
Tc{H) a surface superconducting sheath develops near the sample boundary. The comparison 
of the Tc{H) data obtained on the perforated Al microstructures with that of a reference 
microsquare without antidots confirms the presence of a giant vortex state in the three 
structures in the high magnetic field regime. 

II. EXPERIMENT 

Three different microstructures, shown in Fig. 0, have been studied. A square dot, 
with side a=2.04 fim is taken as a reference sample (a); a square of side a=2.04 fim, with 
four 0.46x0.46 /im^ square antidots (b); and a square, side a=2.14 fim, with two 0.52 x 
0.52 /im^ antidots, placed along a diagonal (c). For the 4-antidot sample (b) the width of 
the superconducting outer stripes is 0.33 fim and the inner stripes are 0.46 fim wide. For 
the 2-antidot sample (c) the outer stripes are 0.35 yum wide, and the non-perforated areas 
are 1.27 /im wide. The dimensions are summarized in Table |. Electrical contacts have been 
attached to the samples using an ultrasonic wire bonding technique on the 150 x 150 /im? 
large contact pads. 

The three samples have been prepared in a single run by thermal evaporation of 99.999 % 
pure Al on a Si02 substrate. The patterns are defined using e-beam lithography on a bilayer 
of PMMA resist previous to the deposition of a 24 nm thick aluminum film. After the evap- 
oration, the liftoff was performed using dichloremethane. The structures were characterized 
by X-ray, SEM and AFM (Fig. p. 

Four-point resistance measurements were performed in a ^He cryostat, using a PAR 124A 
lock-in amplifier. A measuring current of 100 nA r.m.s. with a frequency of 27 Hz was used, 
which is depressing the Tc by only a few millikelvins, in the whole magnetic field range. 

The Tc{H) measurements are done in a continuous run, keeping the sample resistance 
typically at 50 % of the normal state value and sweeping the magnetic field slowly while 
recording the temperature. The magnetic field was applied perpendicular to the structures, 
and a temperature stability better than 0.5 mK was achieved. 
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III. RESULTS 



In Fig. 1^ we present the experimental phase boundary Tc{H) of the three structures. The 
measured Tc{H) values were independent of the direction of the magnetic field scans and 
were reproduced in several measurement rounds. In this paper we will always plot Tc{H) 
in the usual way, i.e. with the Tc-axis pointing from the highest to the lowest temperature. 
Peaks in the Tc{H) plots are then in reality local minima of the critical temperature T^. 

For the reference full square, we observe pseudoperiodic oscillations in Tc{H) super- 
imposed with an almost linear background, where the period of the oscillations slightly 
decreases with increasing field, in agreement with previous studies0S'@ffl. These obser- 
vations are characteristic for the presence of the giant vortex state. For the perforated 
microstructures, two different magnetic field regimes can be distinguished. At high magnetic 
fields, the oscillations in Tc{H) are pseudoperiodic, just as the Tc{H) of the full square. For 
the low field pa.Tt of the phase diagram, distinct features appear (i.e., below ~ 2.5 mT): for 
the 2-antidot sample we observe the same number of peaks compared to the full square, but 
with a considerable shift of the positions of the first peaks. Compared to the full square 
Tc{H), a new series of peaks, positioned symmetrically with respect to /lo-ff ~ 1.4 mT, is 
found for the 4-antidot sample, as can be expected for a 2x2 cell network. 

In what follows, we will investigate in detail the shape of Tc{H) in the two flux regimes 
for the three structures. We will discuss our results in terms of the existing models, within 
the Ginzburg-Landau (GL) theory, developed for mesoscopic structures with a cylindrical 
symmetry (disks, loops) which have been successfully applied earlier to interpret the results 
obtained in mesoscopic square structures0. 



IV. DISCUSSION 

A. The Tc{H) phase boundary of the full microsquare 

The Tc{H) curve measured for the full square structure is very similar to the result 
obtained from a calculationli§ for a mesoscopic disk in the presence of a magnetic field 

(indicated as '-ffcs' in Fig- 0a). In that model the linearized first GL equation is solved with 
the boundary condition for an ideal superconductor/insulator interface: 



(■ 



-ihV - 2eA] ^ 



0, (1) 



which is the condition that no supercurrent can flow perpendicular to the interface. In the 
linear approach, the vector potential A is related to the applied magnetic field H through 
HqH = rot A. In order to obtain the solutions which fulfill the boundary condition (Eq. (|l])), 
one has to solve the equation: 

(L - <l>/$o) M(-n, L + 1, <i>/%) - —^M{-n + 1, L + 2, <l>/<l>o) = (2) 

The function M is the so-called Kummer function of the first kind, n is a real number 
depending on the phase winding number L, which has to be obtained numerically by solving 
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Eq. (H). The flux is defined as $ = hqHtiB^, R being the radius of the disk. The Tc($) is 
obtained via the relation: 



'-'''-■'n^'mir (3) 



T,(0) V ' 2) % 

The upper critical field Hc2 for a bulk superconductor is obtained when substituting = 
in Eq. (^), which gives a linear relation between Hc2 and T. However, for a finite size su- 
perconductor, a third critical field Hcz can be found, because the ground state is obtained 
from solutions of Eq. with n < 0. Superconductivity is concentrated near the sample 
edge (for L > 0), while the "normal" core contains one or several flux quanta L$o- This 
quasi-linear critical field Hcz{T) is the analog of the surface critical field for a semi-infinite 
superconducting slab in contact with vacuum (or insulator), where superconductivity per- 
sists in a surface sheath up to magnetic fields HcsiT) ^1.69 Hc2(T) above the upper critical 
fieldli. 

The series of peaks in the Tc{H) curve correspond to transitions between states with 
different angular momenta L ^ L + 1 of the superconducting order parameter as successive 
flux quanta, $ = L $0) enter the superconductor. A comparison with the experimental result 
for a square structure was made in Ref.0. In a very recent paper by Jadallah et a/.@ the 
Tc{H) phase boundary is studied theoretically and is compared to the experimental Tc{H) 
curve for the full square, described in the present paper. 

Following Ref.cJ, between $=0 and the first peak located at $ = 1.92 $0, the super- 
conducting order parameter \E' has angular momentum L = 0, and the reduced critical 
temperature is quadratic in $: 



T,(0) 2/?2 Uo 

The quasi-linear background at high flux, $/$o ^ 1? follows the asymptotic expression: 

re('i') _ 2e^(o) $ 

T,(0) 71 i?2 



(4) 



(5) 



The parameter rj represents the ratio of the ground state energy (-ffcs) to the lowest bulk 
Landau level (-f^c2) and therefore coincides with the ratio Hc3/Hc2 at a fixed temperature. 
For $/$o ~^ oo, in other words for R oo, the value rj 1.69. Note that substituting 
= 1 in Eq. (|^) gives the equation for Hc2(T) (or Tc2{H)), indicated by the straight hne 
labeled 'H^2' in Fig. |a. 

When matching the position of the experimental and theoretical peaks, we obtain a value 
for the field corresponding to one flux quantum, /io Hq = 0.53 mT. This scaling leads to a very 
good agreement in the position of all the peaks (see the inset of Fig. ||) and strongly supports 
the validity of applying this model to our experiments. From Hq we obtain an effective area of 
3.9 /im^, close to the actual size of the structure, 4.2 /im^. The introduction of this 'effective 
area' is obviously not needed if the Tc{H) is compared with a calculation performed for 
a squareii. Then, from Eq. (^), we find the coherence length, ,^(0)=92 nm (dashed line 
in Fig. Ha), and using, the values for pure Al and the Ginzburg-Landau expressions for 
dirty superconductorsE3, we can estimate the mean free path, i=7 nm, and the penetration 
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depth, A(0)=140 nm. The results for the three structures are summarized in Table |I[ The 
determined value ^(0)=92 nm might be a bit too low, since at low magnetic fields the 
electrical leads attached to the square can give rise to nonlocal effectH. Simultaneously, 
fitting the low field part of the experimental Tc{H) to Eq. (6) of Ref.ii gives an increased 
^(0)=95 nm. 

In contrast to the experimental result presented in Ref.0, which was obtained for a 
substantially larger, but circular dot, the field period can be matched to the theoretical 
predictions in the whole field interval. The distance A<l> between the peaks in Tc($) follows 
the asymptotic limit A$ = $o (l + (2 ?7 $/$o)~^^^) • In our experiment A$ reaches a nearly 
constant value for $/$o > 6, with fioAH ~ 0.60 — 0.65 mT. When a sufficiently high 
magnetic field is applied to the sample, a superconducting edge state is formed, where 
superconductivity only nucleates within a surface layer of thickness wh = ^/^of^V^f^oH. 
The remaining area acts like a normal core of radius RqQ ~ R—wh, and carries L flux quanta 
in its interior. Due to the expanding normal core, the sample can be seen topologically as 
a loop of variable radius. For this reason, the Tc{H) of the dot shows Little-Parks- like 
oscillations, which are, however, nonperiodic. The magnetic period AH decreases, since the 
'effective' radius grows with increasing field. In contrast to the case of the loop, which has 
a parabolic background on Tc{H), the background is quasi-linear, because of the additional 
energy cost (i.e. extra reduction of T^) for depressing superconductivity in the sample core. 
As the applied magnetic field grows this 'giant vortex' core expands until it almost fills the 
entire sample area. According to this expression, at /iq H=5 mT, for example, wh ~ 0.2 //m 
and the effective area of the normal core for the full square structure is ~ 3.3 fim'^. This 
value is in agreement with the observed magnetic period /io AH. 

The amplitude of the experimental oscillations is higher than expected from the theory 
(which was observed also in Ref. 0) (see Fig. I a). We carried out a few Tc{H) measurements, 
where we fixed the electronic feedback circuit at a different resistance value (10 — 90 % of 
the normal state resistance Rn)- When a higher fixed resistance value was chosen (90% of 
Rn), the amplitude of the Tc{H) oscillations was decreased. We should mention here that 
the measured resistance versus temperature curves, as the magnetic field increases, reveal 
a 'bump' above the normal state resistance, together with a significant broadening of the 
temperature interval in which the transition takes place. In Ref.o the bump was shown to 
appear in the 'superconducting state'. Since this disturbance is still small at low fields, we 
believe that the determination of ^(0) from the low field formula (Eq. (^)) is rather accurate 
and therefore we use ^(0)=92 nm further on in this paper. The resistance peaks are caused 
by the electrical leads attached to the square, which have a higher transition temperature. 
For this reason, normal/superconducting boundaries are created at T ?s Tc, which, as in the 
experiment of Park et alB^, can give rise to a resistive transition showing a peak. 



B. The Tc{H) phase boundaries of the perforated microsquares 

Now, we will analyze the phase boundary observed for the 4-antidot structure (Fig. ^o). 
It shows similarities with the full square: there is a quasi-linear background with pseudoperi- 
odic oscillations at high fields and the Tc{H) is parabolic for low magnetic fields (dashed line 
in Fig. 0b). At the same time new features (extra peaks) are clearly seen below ~2.5 mT. 
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Let us first discuss the Tc{H) curve in this low field regime. We will compare it with 
the Tc{H) calculated for a 2 x 2-cell network consisting of one-dimensional strips. Strictly 
speaking, the theory for networks is valid only when the width of the strips forming the 
structure are much smaller than ^(T). For the dimensions of the sample studied here, 
variations of \E' along the strip width can be expected if T is slightly below Tc(0). The de 
Gennes-Alexander (dGA) mo dei, based on the linearized Ginzburg-Landau equations, has 
been used successfullv to explain the phase boundaries obtained in mesoscopic single- and 
multiloop structurescl'Q with narrow superconducting strips. The depression of Tc{H) can be 
expressed as the sum of a topology dependent oscillatory component and a parabolic term 
(dashed line in Fig. Hb) which is due to the finite width w of the strips: 



Tc(0) 3 I $0 J 



(6) 



In Fig. ^ we compare the low magnetic field part of the experimental Tc{H), where a 
parabolic background (Eq. (^) with an averaged (over inner and outer strips) width of 
0.4 fim has been subtracted, with the Tc{H) obtained from the dGA model. In Ref.0 
(Eqs. 20-22) and Ref.i (Eq. 3.12) the functions forming the Tc{H) for a 2 x 2-cell network 
can be found. The field corresponding to one fiux quantum per elementary cell is 2.8 mT, 
leading to an effective total area for the 2 x 2-cell network of 2.96 //m^ (0.74 /im^ per cell, the 
side length of each cell a=0.86 /im). The theoretical Tc{H) reproduce the observed fiux (or 
fiuxoid) states, and have been calculated with ^(0)=92 nm (dashed line in Fig. obtained 
for the full square, and with ^(0)=140 nm (solid line in Fig. estimated from the parabolic 
background (Eq. (|^)). 

It is clear that, for increasing magnetic field, the Tc{H) of the 4-antidot structure can no 
longer be calculated from the dGA model, since ^(T) becomes comparable to the width of 
the strips, giving rise to spatial variation of \l/ perpendicular to the strips. 

Let us point out the differences between the Tc{H) of the present 4-antidot structure and 
the previously measured '2 x 2-antidot cluster' made of Pb / Cu3i. In the present Al structure, 
the network features are not periodically repeated as the magnetic field is increased. Instead, 
above 2.8 mT, the positions of the successive peaks coincide with the peaks observed in the 
Tc{H) of the reference full square. Moreover, the background clearly starts deviating from 
parabolic to quasi-linear. Contrary to this, in the Pb/Cu antidot cluster (see Fig. 1 in 
Ref.i) the peaks related to the network behavior are visible over two periods, i.e. up to 
IIqH ^ 5 mT. For higher fields, no giant vortex state can be deduced for the Pb/Cu sample, 
since the background reduction of Tc stays parabolic and Tc{H) shows pronounced peaks 
instead of cusps. 

The main parameter which determines the Tc{H) is the coherence length ^(T). Since the 
coherence length ^(0) of Al is approximately three times larger than for Pb/Cu, the relative 
Tc reduction 6Tc = 1 — Tc{H) / Tc(0) is almost a factor 10 higher in Al than in an identical 
Pb/Cu sample (see Eqs. (|5)-(0)). Since, for a particular sample geometry, (5Tc/^^(0) = 
1 / ^^(T) should only depend on the magnetic field (at T = Tc{H)), the penetration depth 
A(T) might play an important role. Using A^(T) = A^(0) / ST^, with A(0)=140 nm for Al, 
and A(0)=76 nm for Pb/Cu@'i, we obtain A (T = Tf) ^ 0.6 A (T = T^^/Cuy Q^her 
words, the assumption ^qH = rot A (corresponding to \ ^ w) is fulfilled up to higher 
magnetic fields in the case of Pb/Cu. This is consistent with the fact that, the peaks due to 
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switching of the state of single antidots, is seen up to higher fields in Pb/Cu, but it does not 
explain the very different behavior in the two materials of course. Other possibilities might 
be related to the proximity effect in Pb/Cu, as well as to a different saturation numb ei§ 
ris ^ R/ {2^{T)) in the two materials, although the simple formula for rig was obtained only 
for a single antidot surrounded by a large superconducting area and might not be valid here. 

For the (slightly larger) 2-antidot Al structure (Fig. ^) the interpretation of the low 
field regime of Tc{H) is more difficult. We will return briefly to this point later in the 
paper. At high fields, however, the positions of the peaks in Tc{H) correspond to the same 
pseudoperiodic oscillations as for the full square and 4-antidot structure. 

In Fig. ^ we have replotted the phase diagram in units of $/$o- It is important to 
note that we defined the flux as $ = fioHSeff, with Seff the effective area of the whole 
microsquare. It is close to the exact outer sample area S, and was introduced in order to fit 
the peak positions to the calculated Tc{^) for a circular dot. To avoid confusion, we want 
to mention the different definition of flux in Refs.Bi (antidot cluster), where flux is referred 
to the area available for one single antidot. In the case of a loop0ii, it is natural to define 
the flux as the area enclosed by a contour through the middle of the strips multiplied by the 
magnetic field, in order to ensure a perfectly periodic Tc($), with a period $o- 

Since the 2-antidot structure is a bit larger than the two other structures, a different Hq 
is used to scale the magnetic field. In the inset the positions of the peaks in the experimental 
Tc($) are compared with the theoretical prediction for a mesoscopic superconducting disk. 
The n-th peak corresponds to the transition between the state L = n — 1 and L = n (for 
the 4-antidot sample the peak numbers have been reassigned due to the extra peaks in 
the network regime). At high magnetic fields, there is a quite good agreement of the peak 
positions of the three structures and a good correspondence with the theoretical values found 
for the disk, which is drawn in the inset of Fig. ^ as a solid line. 

How can we understand this striking coincidence of the peak positions at high fields for 
the three structures? For this, we have to look how the superconducting order parameter 
nucleates along a curved superconductor/insulator boundary. Figure ^ shows the calculated 
Tc($) curves for a single circular dot and for an antidot (see also Ref.El) in an infinite film, 
both of radius R. The latter has been calculated in a similar way as the Tc($) of the dot. 
For a single antidot, the boundary condition (Eq. (|l])) translates into: 

(L - $/<l>o) Ui-n, L + 1, $/$o) + 2n-^ U{-n + 1, L + 2, $/<l>o) = , (7) 

where the function U is the Kummer function of the second kind, diverging at the origin, 
i.e. at the center of the antidot. The numeric values for n have to be inserted into Eq. (H), 
to obtain the Tc($). 

In Figure the respective enhancement factors t] (corresponding to Hc3/Hc2) are shown. 
For both the dot and the antidot the value 77=1.69 is approached as the curvature radius 
R goes to infinity (or for a fixed R, a.s H oo). Since the dot has a larger rj than the 
antidot, corresponding to a higher Hc3{T), the superconducting order parameter is expected 
to grow initially at the outer sample boundary, as the temperature drops below T^. At 
slightly lower temperatures surface superconductivity should as well nucleate around the 
antidots. In the mean time, however, the order parameter has reached already a finite value 
over the whole width of the strips. In the complete temperature (or flux) interval of our 
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measurements rj < 1.5 for the antidots and f] > 1.8 for the dot (when scahng the radii 
to the actual sample dimensions). The resistively measured phase transitions, probably 
because of this substantially different H^s for a dot and an antidot, only show peaks related 
to the switching of the angular momentum L, associated with a closed contour along the 
outer sample boundary. At the Tc{H) boundary, in the high magnetic field regime, there is 
no such closed superconducting path around each single antidot, and therefore the fluxoid 
quantization condition does not need to be fulfilled for a closed contour encircling each single 
antidot. 

Although a more detailed analysis has to be carried out, since the sample boundaries in 
our experiments have sharp corners, the interpretation given above is expected still to be 
valid. A detailed analysis of a square loop geometry performed by Fomin et a/.ii has shown 
that the superconducting order parameter preferentially nucleates near the sharp corners of 
the structure. In another paper, the same authors discuss the enhancement of the surface 
critical field H^s above the bulk upper critical field Hc2 in a semiplane, which is bent over 
a certain angle a (superconducting wedge )!!. The magnetic field is parallel to the wedge 
edge. An enhanced rj value is found for angles a < tt , which can be as high a.s 1] = 3.79 
for a = 0.44 TT. For angles a > it, the surface critical field is not enhanced above r] = 1.69. 
Note that the value t] = 3.67, obtained for a = 7i/2 differs from the calculation in Ref.@, 
where for a square domain in the limit $/$o cx) the factor 1.8 only. The discrepancy 
between these two results still has to be clarified. 

For the 2-antidot structure, further quantitative calculations are needed to describe the 
low field behavior. Since there are no extra peaks present in Tc($), compared to the full 
square, we believe that, here as well, a surface superconducting sheath develops along the 
outer sample boundary. At the lowest fields, the sheath width wh is still larger than the 
width w of the strips, and therefore the position of the first peaks (mainly the second) in 
Tc($) are different from the full square. 

In Fig. ^ the superconducting order parameter profiles are shown for a disk, calculated 
at different points $/$o = 1,3, 17 on the Tc($) curve. For $/$o = 1 the ground state 
corresponds to L = 0, and the order parameter is only weakly modulated since there are 
no flux lines threading the sample. As we move to higher $/$o superconductivity becomes 
more and more concentrated near the sample boundary. The presence of the antidots in 
the perforated samples produces different profiles for the superconducting order parameter, 
since the \E' has to fulfill the boundary conditions (Eq. (|l])) also at the antidot boundaries. 
A two-dimensional GL calculation would be required to obtain the proper order parameter 
distribution here. 

The inset of Fig. |^ shows the width wh of the edge state as a function of $/$o- For 
our samples wh becomes equal to the width of the strips at ~ (3 — 4) $/$o- For fluxes 
above this value the presence of the antidots will not influence the position of the peaks 
in Tc($). Therefore, at high fields, when wh is smaller than the width of the outer strips 
in our structures, the order parameter is restricted to the outer border of the sample, and 
therefore it is impossible to have supercurrents around a single antidot. 

The background depression of Tc is different for the three structures studied (Fig. |). The 
larger the perforated area (in other words the smaller the area exposed to the perpendicular 
magnetic field), the less Tc($) is pushed to lower temperatures. Another clear example of 
a similar behavior is given in Ref.0, where the Tf.{H) of the (square) dot is shown to be 
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lower than the Tc{H) of a the loop, when exposed to a perpendicular magnetic field. This 
general rule applies, for instance, also to simple strips for which Tc{H) is suppressed more 
when the width w increases, which is described by Eq. (^). In the dot of Ref.§, the giant 
vortex state was shown to develop. For the loop studied in that paperS, however, the used 
magnetic fields were too low to induce the crossover to a giant vortex state, which contrasts 
the observations from the present paper, in the high magnetic field regime. 

The appearance of the giant vortex state in the high field regime is the most plausible 
explanation at the moment. We can not exclude, however, that another scenario is also 
possible. Namely, a nearly fiat non-zero distribution of l^'l in the sample interior could 
coexist with an enhanced at the external sample boundary, although we believe that 
such a situation would give rise to peaks in Tc($) each time an antidot changes its quantum 
state. In any case, the final description of the specific shape of the superconducting order 
parameter at Tc($) requires a numerical two-dimensional calculation of \l/ for the perforated 
topologies, where the boundary conditions are fulfilled both at the outer and at antidot 
superconducting/ insulator interfaces. 

In summary, we have presented the experimental superconducting/normal phase bound- 
aries Tc{H) of a mesoscopic full square and two perforated mesoscopic aluminum squares. 
The fiux interval was divided in two regimes by comparing the results with the behavior of 
a full square microstructure: for low magnetic fields the 4-antidot structure behaves like a 
network consisting of quasi-one-dimensional strips, giving rise to extra peaks in Tc{H) in 
comparison to the full square. In the 2-antidot structure the peak positions are only shifted 
compared to the full square. As soon as each antidot contains one fiux quantum, the gi- 
ant vortex develops, resulting in pseudoperiodic oscillations in the Tc{H) and a quasi-linear 
background on Tc{H) at high magnetic fields. In this regime, the peak positions coincide 
for all three structures studied when the phase boundaries are plotted in fiux quanta units 
(where fiux is referred to the total sample area). For high magnetic fields, the presence of 
the antidots apparently does not change the phase winding number L for a closed contour 
around the outer perimeter of the whole square. Since the enhancement factor Hc3/Hc2 
is the highest at the outer sample boundary, superconductivity nucleates initially near the 
outside sample edges, resulting in a giant vortex state. 

Note added invroof: The discrepancy between the result obtained by the authors of 
Ref.ii and of Ref.ii has been clarified in an erratumB. 
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FIGURES 

FIG. 1. AFM images of the three structures: (a) full, (b) 4-antidot, and (c) 2-antidot mi- 
crosquares. 

FIG. 2. Tc{H) phase boundaries for the (a) full, (b) 4-antidot, and (c) 2-antidot microsquares. 
For the full microsquare (a) we also present the calculated Hcs for an equivalent disk (Eq. (|3|)). 
The straight solid line is the calculated upper critical field Hc2 (substitute n = in Eq. (^). The 
parabola (dashed) indicates the low field Tc{H) behavior (Eq. (Q)). The dashed line in (b) gives 
the parabolic background depression of (Eq. (^)) for structures made of quasi-one-dimensional 
strips. 

FIG. 3. Tc{^) phase boundaries in reduced units of critical temperature and flux. For $/$o > 5 
the peaks in Tc(<I>) appear at the same flux $/$o in ^ the structures. The inset shows a comparison 
of the peak positions in Tc{^) with the predictions for a mesoscopic disk. The scaling of the 
theoretical values (in flux units) gives the effective area of the microsquare. 

FIG. 4. Low field part (single period) of the experimental phase boundary Tc{H) of the 4-an- 
tidot sample (where a parabolic background (Eq. (^) has been subtracted) compared with the 
Tc(-fr) calculated from the dGA model for a 2 x 2 cell network made of one-dimensional strips (see 
Refs.0jl), and using ^(0)=92 nm (dashed line), obtained for the full square, and ^(0)=140 nm 
(solid curve), estimated from the parabolic background (Eq. (P)) and an average strip width of 0.4 

FIG. 5. (a) Calculated phase boundaries (i.e., the Hcz{T) curves) for a circular antidot (Eq. (|^) 
and a dot (Eq. (^)) in normalized units of temperature and magnetic flux. Superconductivity will 
always nucleates initially near the dot /insulator boundary (the dot has the highest i/cs)- (b) Plot 
of the enhancement factor r/ (corresponding to Hc3/Hc2) for the same structures. 

FIG. 6. Calculation of the modulus of the superconducting order parameter |^'| for a disk, 
where R is the disk radius, and r the distance measured from the sample center to the sample 
boundary. The different curves are calculated for several values for $/$o ^-t the phase boundary 
Hc2,{T). The inset shows the normalized width of the superconducting edge state as a function of 
^>/$o- 
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TABLES 



Low H 


area (/um^) 


TM (K) 


fioHo (mT) 


eff. area (/^m^) 


C(O) (nm) 


i (nm) 


A(0) (nm) 


full square 


4.16 


1.361 


0.53 


3.9 


92 


7 


140 


4-antidot 


4.16 


1.353 


0.53 


3.9 


92 


7 


140 


2-antidot 


4.58 


1.369 


0.57 


4.15 


92 


7 


140 



TABLE L Values of several characteristic magnitudes for the three structures. The field cor- 
responding to one flux quantum, ^iqHq, and the effective areas are obtained by matching of the 
position of high field peaks to the theoretical location for a mesoscopic disk. The characteristic su- 
perconducting lengths for the full square are obtained from the low field behavior of Tc{^) (Eq. (Q)) 
and we assume the same values for the other structures. 



14 



Fig.2 




1 2 3 4 5 6 7 



Hj,H(mT) 



Fig. 5 




